Back

Genes & Immunity

Springer Science and Business Media LLC

Preprints posted in the last 90 days, ranked by how well they match Genes & Immunity's content profile, based on 11 papers previously published here. The average preprint has a 0.01% match score for this journal, so anything above that is already an above-average fit.

1
Large-scale meta-analysis of over one million individuals reveals the genetic architecture of 127 complex traits in East Asian populations

Jo, J.; Khor, S.-S.; Chu, S.-K.; Ji, Y.; Ueno, K.; Ono, A.; Chen, C.-W.; Do, A.; Han, H.; Kawai, Y.; Kim, N.-E.; Chen, C.-h.; Tokunaga, K.; Won, S.; Yang, H.-C.

2026-06-23 genetics 10.64898/2026.06.18.730290 medRxiv
Top 0.1%
2.4%
Show abstract

Genome-wide association studies (GWASs) have disproportionately focused on European (EUR) populations, limiting the characterization of genetic architecture in other ancestries. To address this imbalance, we integrated large-scale biobanks from Japan, Korea, Taiwan, and China to perform the largest phenome-wide meta-analysis to date in East Asian (EAS) populations, encompassing over one million individuals across 127 complex traits. We identified 8,010 previously unreported associations and observed substantial genetic sharing across EAS subpopulations, while also detecting cohort-specific heterogeneity within the broader EAS context. Transethnic analyses revealed moderate genetic correlations between EAS and EUR populations, indicating both shared and ancestry-specific components of disease risk. Pleiotropy analyses highlighted prominent signals within the HLA region, supported by protein-protein interaction connectivity and immune-related pathway enrichment. Decomposition of genome-wide association matrices further uncovered structured cross-trait architectures, revealing a predominantly shared polygenic backbone driven by metabolic, biochemical, and anthropometric traits, together with two discrete latent components enriched for immune-related processes. Together, our findings refine the genetic architecture of complex traits in East Asian populations at unprecedented scale and clarify the balance between shared and population-specific determinants of human diseases.

2
Genome-wide meQTL mapping in cattle blood reveals cis and trans regulation of DNA methylation

Fouere, C.; Costes, V.; Besnard, F.; Le Danvic, C.; Patry, C.; Fritz, S.; Boussaha, M.; Jouin, M.; Boichard, D.; Kiefer, H.; Costa Monteiro Moreira, G.; Sanchez, M.-P.

2026-07-08 genetics 10.64898/2026.07.07.736355 medRxiv
Top 0.1%
2.4%
Show abstract

Background Complex traits are influenced by numerous variants, most of which have regulatory effects on gene expression that can be mediated by DNA methylation. Molecular QTL mapping is an approach that aims to dissect these effects. However, obtaining molecular phenotypes on a large scale is challenging, particularly in livestock species. In cattle, an epigenotyping array called EpiChip has recently been developed in the European RUMIGEN project. The EpiChip, which contains 43,317 CpG sites distributed all over the bovine genome, enables large-scale measurement of DNA methylation. This study aims to characterize the genetic determinism of blood DNA methylation in cows by estimating heritability and mapping cis- and trans-methylation QTLs (meQTLs). Results Whole blood samples from 4,457 genotyped Holstein cows were epigenotyped. Across all CpG sites, the heritability estimates averaged 24.6%. The local meQTL mapping at sequence-level for variable CpG sites (SD > 2.5%; n = 28,806) detected cis-meQTLs for 80.1% of the CpG sites, with sentinel SNPs located close to their associated CpGs. A two-step analysis was also conducted to identify long-range associations, with a particular focus on trans-meQTL hotspots. First, we identified CpG-SNP trans-associations using medium-density genotypes (50k SNPs) that revealed 31,846 SNPs with significant effects on 1 to 530 trans-CpG sites. Then, regions associated with at least 34 independent trans-CpGs were retained defining 31 hotpots. For each hotspot, a local sequence-level GWAS was conducted using the first principal component derived from the associated trans-CpGs. Out of the 31 detected hotspots, three were located close to transcription factor genes (RUNX1, NFIC and FOXA3) for which the associated trans-CpGs were enriched for the corresponding binding motif. Two other hotspots were located within KDM5A and KDM5B, and their corresponding trans-CpGs were strongly overrepresented in H3K4me3 narrow peaks in blood as well as in other tissues. Conclusions By identifying functional candidate genes associated with blood DNA methylation in cattle, these findings provide new insights into the regulatory architecture of DNA methylation in mammals, highlighting the value of large-scale molecular data from livestock populations.

3
A systematic analysis of IBD GWAS loci identifies most probable causal genes impacting intestinal epithelial functions

Hebert-Milette, I.; Mercier, V.; Paquette, J.; Boucher, G.; Levesque, C.; Goyette, P.; Rioux, J. D.

2026-08-02 genetic and genomic medicine 10.64898/2026.07.30.26359349 medRxiv
Top 0.1%
1.7%
Show abstract

Background Genome-wide association studies have identified >200 loci associated with IBD, yet the causal gene for most remains unknown. As multiple epithelial functions have been linked with susceptibility to IBD, there is a need to prioritize candidate causal genes for functional studies in this cellular context. Methods Using a standardized definition of regions implicated by index SNPs from three GWAS studies, we categorized regions as containing: (1) a known casual gene, (2) a single gene or (3) multiple genes. We then developed an IBD Priority Score to rank genes based on genetic, genomic and functional data. We next developed and applied an Epithelial Priority Score, based on expression patterns and quantitative traits, to prioritize genes for functional validation in epithelial models. Two candidate genes identified through this approach were tested for their impact on viral response pathways in HT-29 cells. Results The IBD Priority Score prioritized a single gene in 71 of the 104 regions containing multiple genes. The Epithelial Priority Score identified 31 epithelial candidates. Functional studies demonstrated that IRF6 enhanced, whereas IRF8 suppressed, antiviral responses in intestinal epithelial cells stimulated with Poly(I:C). Conclusions Combining multiple genetic, genomic, and functional data is a useful approach for prioritizing the most likely causal gene within IBD GWAS loci, and for prioritizing functional validation studies in epithelial cells and tissues. Moreover, we provide functional evidence for two IBD genes playing a role in the regulation of anti-viral responses in intestinal epithelial cells.

4
Maintenance DNA methylation is necessary for age-related alterations in regulatory T cell transcriptional and DNA methylation signatures

Gurkan, J. K.; Liu, Q.; Reyes Flores, C. P.; Helmin, K. A.; Ryan, D. H.; Joudi, A. M.; Ulrich, B. J.; Abdala-Valencia, H.; Steinert, E. M.; Singer, B. D.

2026-06-29 immunology 10.64898/2026.06.24.733005 medRxiv
Top 0.1%
1.7%
Show abstract

CD4+FOXP3+ regulatory T (Treg) cells maintain self-tolerance, restrain immune responses during inflammatory stimuli, and promote tissue function and repair. Treg cell lineage identity, stability, and function depend on specific DNA methylation patterns maintained by the epigenetic regulator, UHRF1. Aging disrupts DNA methylation patterns necessary for Treg cell-mediated lung repair in a cell-autonomous manner. Nevertheless, whether maintenance DNA methylation is necessary for age-related Treg cell transcriptional and methylation programs is unknown. Here, we performed transcriptional and DNA methylation profiling on young and old Treg cells isolated from mice with chimeric Treg cell-specific loss of UHRF1. We observed cell-autonomous, age-related alterations in transcriptional and DNA methylation signatures that were dependent on UHRF1. We conclude that maintenance DNA methylation is required for age-related alterations in Treg cell transcriptional and DNA methylation signatures.

5
Epigenetic dysregulation of Th2 cytokine genes in MuSK myasthenia gravis and its modulation by immunosuppressive therapy

Elmas, C.; Stoccoro, A.; Lari, M.; Salehi, F.; Iovino, V.; Cepele, A.; Huber, J.; Faber, F.; Wolfsgruber, M.; Keritam, O.; Weng, R.; Steinmaurer, A.; Koenig, T.; Guida, M.; Cetin, H.; Zimprich, F.; Hoeftberger, R.; Maestri Tassoni, M.; Coppede, F.; Koneczny, I.

2026-08-11 immunology 10.64898/2026.08.05.742975 medRxiv
Top 0.1%
1.5%
Show abstract

Background and objectivesMyasthenia gravis associated with antibodies against muscle-specific kinase (MuSK-MG) is a well-characterized IgG4-autoimmune disease, however, the mechanisms driving IgG4 predominance remain poorly understood. This study investigated whether promoter DNA methylation of cytokine genes involved in IgG4 class switching is associated with this immune response. MethodsPeripheral blood mononuclear cells were isolated from MuSK-MG patients (n=36), acetylcholine receptor myasthenia gravis (AChR-MG) patients as disease controls (n=7), and sex-matched healthy controls (n=12). Promoter DNA methylation of IL4, IL10, and IL13 was assessed by methylation-sensitive high-resolution melting and relative cytokine mRNA expression by qPCR. Associations with clinical variables, and antibody levels were subsequently evaluated. ResultsMuSK-MG patients showed lower median IL13 promoter methylation compared with healthy controls (p = 0.004). Median IL4 promoter methylation was also reduced in MuSK-MG compared with healthy controls (p < 0.001) and AChR-MG disease controls (p < 0.001), whereas no differences were observed for IL10 promoter methylation. Relative mRNA expression of IL4 (p = 0.0005), IL10 (p = 0.0462), and IL13 (p = 0.0002) was increased in MuSK-MG compared with AChR-MG. Compared with healthy controls, only IL4 expression remained significantly increased (p < 0.0001). Promoter methylation was inversely correlated with relative mRNA expression for IL4 (p < 0.0001), while IL13 showed a similar but non-significant trend (p = 0.054), no association was observed for IL10. Multivariable analysis demonstrated that treatment at sampling was independently associated with lower IL10 and IL13 promoter methylation, whereas no associations were observed with age, sex, disease phase, or disease duration. Promoter methylation did not correlate with total serum IgG4 or anti-MuSK IgG4 levels. DiscussionMuSK-MG is associated with selective hypomethylation of IL4 and IL13 promoters accompanied by increased cytokine gene expression, while IL10 promoter methylation remains unchanged. The association between treatment and IL10 and IL13 promoter methylation suggests that immunosuppressive therapy may influence epigenetic regulation in MuSK-MG. Together, these findings support a role for epigenetic dysregulation of Th2-associated cytokines in the immunological environment associated with IgG4 subclass switch. To our knowledge, this is the first study investigating IL4, IL10, and IL13 promoter DNA methylation in MuSK-MG.

6
Placental DNA methylation captures shared and trait-specific genetic susceptibility across complex health conditions

Cilleros Portet, A.; Gonzalez-Moro, I.; Saddiki, H.; Mari, S.; Everson, T.; Hernangomez-Laderas, A.; Broseus, L.; Tost, J.; Cosin-Tomas, M.; Groleau, M.; Czamara, D.; Lozano, M.; Hao, K.; Tuhkanen, J.; Fallin, M. D.; Schmidt, R. J.; Breeze, C. E.; Deleuze, J.-F.; Aguilar-Lacasana, S.; Jacques, P.-E.; Hytti, S.; Irizar, A.; Lahti-Pulkkinen, M.; Bakulski, K. M.; Dou, J.; Lahti, J.; Vrijheid, M.; Raikkonen, K.; Hivert, M.-F.; Sunyer, J.; Heude, B.; lepeule, j.; London, S. J.; Chen, J.; Bustamante, M.; Marsit, C.; Bilbao Catala, J. R.; Lesseur, C.; Fernandez-Jimenez, N.

2026-07-09 genetic and genomic medicine 10.64898/2026.07.07.26357471 medRxiv
Top 0.1%
1.5%
Show abstract

The Developmental Origins of Health and Disease (DOHaD) hypothesis proposes that the perinatal environment shapes susceptibility to complex traits across life [1]. The placenta, a transient organ mediating maternal-fetal exchange, plays a central role in this process and has emerged as a key molecular archive in utero [2-4]. Placental DNA methylation (DNAm) is a unique mediator between prenatal exposures, fetal genetics and later-life outcomes [5-9]. DNAm quantitative trait loci (mQTL) have helped disentangling causal mechanisms underlying GWAS loci for complex diseases [10-15]. Despite growing evidence that placental genomic regulation has broad and profound effects on the developmental programming of early- and later-life health outcomes [17], existing placental studies remain limited in scale and largely focused on growth- and neuro-related traits [12-16]. Here, we construct a high-resolution placental mQTL resource and systematically investigate how placental DNAm relates to early- and later-life traits, and to shared vulnerability and complex interactions among them.

7
Long-range regulatory target prediction reveals shared genetic background across ulcerative colitis, Crohn's disease, primary sclerosing cholangitis and ankylosing spondylitis

Dulcic, D.; Mandic, K.; Hrsak, D.; Baresic, A.

2026-07-03 genomics 10.64898/2026.06.29.735270 medRxiv
Top 0.1%
1.3%
Show abstract

Common variants detected by the genome-wide association studies (GWAS) create a wealth of knowledge on genetic component of individual traits and diseases. Elucidating the molecular mechanism behind the vast majority of these variants that are found to be non-coding remains a largely unsolved task, especially when distal and pleiotropic interactions between regulatory elements where these variants occur and gene promoters are taken into account. Focusing on four diseases with immune-mediated mechanisms namely ulcerative colitis, Crohn's disease, primary sclerosing cholangitis and ankylosing spondylitis, we demonstrate the utility of the targPred tool, providing prediction of genes targeted by the regulatory variants. We demonstrate that taking into account evolutionary and comparative genomic data, previously unobserved mechanistic trends (the platelet, vascular and sterol clusters) can be detected in terms of implicated genes targeted by the regulatory elements containing common variants, shared between all four diseases, as well as specific trends for subsets of diseases, e.g. two IBD phenotypes. We also elucidate a clinically-relevant target COG6 shared between IBD and PSC, as well as a whole range of other target genes missed by the conventional SNP-to-gene assignments methods.

8
3D genome remodelling underlies PU.1-dependent granulocyte maturation

Jia, W.; Johanson, T. M.; Dakic, A.; Garnham, A.; Smyth, G. K.; Nutt, S. L.; Allan, R. S.; Coughlan, H. D.

2026-07-21 immunology 10.64898/2026.07.21.739714 medRxiv
Top 0.1%
1.3%
Show abstract

PU.1 is an ETS-domain transcription factor that has critical roles in many aspects of hematopoiesis and immune cell fate and function. In addition, aberrant PU.1 expression has been implicated in the development of acute myeloid leukemia (AML). Loss of PU.1 during adult murine hematopoiesis results in the expansion of immature granulocytes suggesting that PU.1 plays an important role in granulocyte maturation. To understand the molecular underpinnings of this process, we combined gene expression, transcription factor binding and 3D genome analyses with conditional deletion of PU.1 in vivo. We find that in contrast to normal granulocytes, PU.1-deficient cells possessed a transcriptome of immature granulocytes, in line with their cellular phenotype. Furthermore PU.1-deficient granulocytes display altered 3D genome architecture with a significant loss of interactivity in regions bound by PU.1 in control cells. Overall, this study implicates PU.1 as a key regulator of granulocyte maturation and lineage commitment through control of transcriptional programs and 3D chromatin architecture.

9
Genome-wide association study of susceptibility to pneumococcal carriage amongst children

Kandasamy, R.; Gurung, M.; Shrestha, S.; Bibi, S.; Thorson, S.; Carter, M.; O'Connor, D.; Murdoch, D. R.; Kelly, D. F.; Shrestha, S.; Levin, M.; Pollard, A. J.

2026-07-16 genetic and genomic medicine 10.64898/2026.07.13.26356474 medRxiv
Top 0.1%
1.1%
Show abstract

Background Pneumococcal disease is a leading cause of paediatric pneumonia and meningitis. Pneumococcal colonisation is the fundamental step to pneumococcal disease causation. We aimed to identify genetic loci associated with pneumococcal colonisation amongst children. Methods We conducted a genome-wide association study on 2111 Nepalese children, comprising 1346 cases carrying pneumococcus and 765 controls. We tested 8.1 million imputed variants using logistic regression and ten principal components as covariates. Fine mapping and functional evidence were used to identify suspected causal variants and related genes of interest. Findings A cluster of 22 variants of genome-wide significance (p<5x10-8) were identified on chromosome 12q21.31, eight of which were within PPFIA2. Fine mapping of this region identified 5 variants within 0.1 Mb of the 5-prime region of PPFIA2 all of which are significant eQTLs for PPFIA2. We further describe three loci (10q23.31, 12q23.1, and 20p11.21) which had variants with highly suggestive associations (p<5x10-7)with pneumococcal carriage. Interpretation Our study demonstrate human susceptibility to pneumococcal carriage to be polygenic with genetic variations which regulate PPFIA2 expression playing a key role in the ability for pneumococcus to colonise children. Targeting these genetic factors and the associated pathways are a means for preventing pneumococcal disease. Funding This study was supported by funding from Gavi - the vaccine alliance, the European Unions Horizon 2020 research and innovation program under grant agreement number 668303 (PERFORM), and a Robert Austrian Research Award.

10
Akkermansia muciniphila-derived LPS links gut dysbiosis to pathogenic miR-21 signaling in experimental autoimmune encephalomyelitis.

Mallahalli, M. S.; Hohjoh, H.; Takewaki, D.; Kimura, K.; Oki, S.; Mori, H.; Hosomi, K.; Kunisawa, J.; Toyoda, A.; Sato, W.; Yamamura, T.

2026-07-10 immunology 10.64898/2026.07.06.736880 medRxiv
Top 0.2%
1.0%
Show abstract

Multiple sclerosis (MS) is a chronic T cell-mediated autoimmune disease characterized by blood-brain barrier (BBB) disruption, neuroinflammation, and demyelination of the central nervous system (CNS). Emerging evidence links gut microbiota to disease pathogenesis, but the microbial factors that regulate pathogenic microRNA (miRNA) programs are largely unknown. Here, using experimental autoimmune encephalomyelitis (EAE, a MS mouse model), we investigated whether gut microbiota exacerbate EAE pathogenesis by modulating host miRNA expression. Antibiotic-induced depletion of the gut microbiota markedly attenuated EAE scores and reduced circulating inflammatory miRNAs, with miR-21 emerging as the dominant pathogenic candidate. Functional inhibition of miR-21 significantly ameliorated disease severity and reduced CNS T-cell infiltration. Mechanistically, miR-21 enhanced IL-17 and GM-CSF production by CD4 T cells and promoted immune-cell entry into the CNS through endothelial activation and blood-brain barrier dysfunction. We identified a transient expansion of Akkermansia muciniphila during the prodromal phase of EAE that positively correlated with circulating miR-21 levels. Colonization of antibiotic-treated mice with A. muciniphila exacerbated EAE and increased serum miR-21, whereas monocolonization of germ-free mice was insufficient to induce systemic miR-21, indicating a requirement for an inflammatory host environment. Further analyses revealed that atypical lipopolysaccharides (LPS) derived from A. muciniphila induce epithelial miR-21 production through coordinated TLR2/TLR4 signaling. Circulating miR-21 subsequently promoted endothelial dysfunction through the TIMP3-ADAM17 pathway, facilitating pathogenic T-cell migration into the CNS. Importantly, circulating miR-21 was also elevated in patients with MS. Collectively, these findings identify a previously unrecognized A. muciniphila-LPS-miR-21 axis linking gut dysbiosis to neuroinflammation and suggest that host-derived miRNAs function as systemic mediators through which microbial signals influence CNS autoimmunity.

11
Liver microbiome composition associates with histological severity and PNPLA3 genotype in metabolic dysfunction-associated steatotic liver disease

Mascardi, M. F.; Taussig, R.; Signoretta, I. P.; Suarez, B.; Marciano, S.; Casciato, P.; Narvaez, A.; Haddad, L.; Gadano, A.; Penas-Steinhardt, A.; Bustamante, J. P.; Trinks, J.

2026-07-09 molecular biology 10.64898/2026.06.30.735597 medRxiv
Top 0.2%
1.0%
Show abstract

BACKGROUNDMetabolic dysfunction-associated steatotic liver disease (MASLD) is a systemic immunometabolic disorder rapidly increasing worldwide, affecting nearly 38% of adults. Gut dysbiosis and host genetic factors, such as PNPLA3 I148M variant, modulate disease development and progression. Through the gut-liver axis, increased intestinal permeability enables microbial translocation to the liver, promoting inflammation and metabolic disruption. However, the composition and functional potential of the hepatic microbiome remain poorly characterized. Understanding its relationship with histological injury and genetic susceptibility may provide novel mechanistic insights. We hypothesized that the hepatic microbiome composition and function are associated with histological severity and PNPLA3 genotype in this disease. AIMTo characterize the hepatic microbiome and assess its association with histological severity and PNPLA3 genotype. METHODSThis cross-sectional observational study included 30 patients with MASLD from a tertiary care hospital. Liver tissue underwent shotgun metagenomic sequencing. Histological severity was assessed using the NAFLD Activity Score (NAS). PNPLA3 genotype was determined by PCR. Differential abundance and functional enrichment analyses were performed using MaAsLin2. Somatic variants were identified using Mutect2. Correlation networks were constructed using Spearmans correlation coefficients. RESULTSPatients with advanced histological injury (NAS [&ge;]5) and PNPLA3 I148M carriers showed a trend toward higher somatic mutational load and a markedly reduced microbial abundance. Analyses revealed broad compositional shifts across bacterial, fungal, viral, and eukaryotic taxa, affecting both commensal and context-dependent pathobiont lineages. Pseudomonas species were enriched, whereas Siphoviridae phages were depleted in advanced disease and PNPLA3 I148M carriers. Functional analysis revealed enrichment of pathways related to nutrient transport and metabolic stress adaptation, while TonB-associated functions were enriched in advanced liver injury but depleted in PNPLA3 I148M carriers. Network analysis identified Sphingomonas leidyi as a keystone node associated with hexosamine metabolism. Salmonella enterica abundance positively correlated with somatic variant burden, suggesting a link between microbial signatures and genomic instability. Histological progression and the risk PNPLA3 genotype were accompanied by marked topological simplification, reflecting less resilient community structures. CONCLUSIONSThe hepatic microbiome in MASLD is a low-biomass, polymicrobial ecosystem shaped by the host genetic background. Its functional activity, taxonomic composition and system architecture bidirectionally relate to liver DNA instability and the severity of histological damage. Core tipThis study characterizes the multi-kingdom hepatic microbiome in MASLD using FFPE-derived metagenomics. We demonstrate that microbial abundance-including bacteria, fungi, protozoa, and viruses- significantly decreases with increased histological severity and the PNPLA3 risk genotype. Rather than global diversity shifts, results showed that disease progression could be linked to specific functional adaptations and simplified microbial network connectivity. In addition, we described associations between specific taxa and somatic mutational burden, suggesting a link between microbial signals and genomic instability. These findings indicate that changes in the liver microbiome as a whole, rather than specific taxonomic modifications, influence MASLD pathophysiology.

12
Single-cell analysis of an adult IBD INCEPTION cohort reveals Galectin-linked disease mechanisms

Leipner, M.; Rimmer, P.; Tull, S.; Paun, A.; Sandrin, V.; Begum, J.; Mansour, A. A.; Saviano, A.; Sharma, N.; Cheesbrough, J.; Maione, F.; Trenkle, P.; Klein, A.; Danilin, S.; Iqbal, T. H.; Iqbal, A. J.; Regan-Komito, D.

2026-07-03 immunology 10.64898/2026.06.30.735473 medRxiv
Top 0.2%
1.0%
Show abstract

Background and Aims: The molecular pathogenesis of Inflammatory Bowel Disease (IBD) remains unclear. We aimed to establish a high-resolution immune landscape of treatment-naive IBD to identify central drivers of disease onset and early pathogenic signalling. Methods: We generated a single-cell atlas using intestinal biopsies from a large adult inception cohort of 137 individuals, including treatment-naive Crohn's disease (CD), ulcerative colitis (UC), and symptomatic non-IBD controls. We integrated scRNA-seq (1 million cells) with co-varying neighbourhood analysis (CNA) and unbiased tensor decomposition of cell-cell communication (CCC) networks. Findings were validated in vitro macrophage stimulation model and using serum from patients. Results: The inception cohort exhibited significantly more homogenous compartmental diversity compared to benchmark reference studies (p < 0.001). Inflammation in both CD and UC was characterized by a marked expansion of inflammatory monocytes. Unbiased CCC analysis identified a dominant disease-specific signalling module centred on the Galectin family (LGALS1 and LGALS9). Galectin-9 expression was specifically enriched in inflammatory monocytes, which exhibited distinct. transcriptional programs linked to antigen presentation and microbial sensing. In vitro, Galectin-9 acted as a potent stimulus, driving macrophages toward a pro-inflammatory phenotype. Clinically, serum Galectin-9 levels were significantly elevated in IBD patients and correlated with systemic inflammatory markers and treatment response. Conclusions: Our data identify a galectin-monocyte signalling axis as a unifying inflammatory hallmark of early IBD. Galectin-9 serves as both a functional driver of mucosal inflammation and a dynamic biomarker, offering new opportunities for therapeutic targeting and disease monitoring from diagnosis. Keywords: Inflammatory Bowel Disease; Crohn's Disease; Ulcerative Colitis; Single-cell RNA sequencing; Galectin-9; Inflammatory monocytes.

13
An integrated human immunoglobulin germline resource linking allele diversity to expressed repertoire structure

Peres, A.; Jana, U.; Rodriguez, O.; Vanwinkle, Z. M.; Engelbrecht, E.; Gibson, W. S.; Shields, K.; Croslin, B.; Schultze, S.; Bharadwaj, C.; Murray, C.; Lees, W. D.; Smith, M. L.; Watson, C. T.; Yaari, G.

2026-06-09 immunology 10.64898/2026.06.05.730236 medRxiv
Top 0.2%
1.0%
Show abstract

Human immunoglobulin (IG) loci are highly polymorphic, yet existing germline resources remain noisy and incomplete, limiting our ability to link inherited variation to antibody repertoires. Here, we integrate high-fidelity long-read genomic sequencing with matched adaptive immune receptor repertoire sequencing (AIRR-seq) to construct HUSA, a population-scale, evidence-resolved germline resource. Using a conservative allele inference framework, HUSA expands current references more than three-fold, identifying over 1300 alleles while preserving allele-level evidence provenance across genomic and repertoire data. By linking genotype and expressed repertoires within individuals, we show that coding-region similarity predicts the structure of adjacent recombination signal sequences and leader regions, revealing that IG alleles are organized as linked cis-regulatory units associated with differences in recombination context and allele usage. These results define key germline constraints shaping repertoire formation and establish a robust, genotype-aware foundation for the analysis of immune receptor repertoires.

14
Bulk and single-cell transcriptomics reveal elevated endogenous retrovirus expression linked to immunopathology in severe COVID-19.

Wang, B.; Deckers, T.; Liu, E.; Tokuyama, M.

2026-07-09 immunology 10.64898/2026.07.08.737320 medRxiv
Top 0.2%
1.0%
Show abstract

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection triggers expression of endogenous retroviruses (ERVs), but whether this persists during coronavirus disease 2019 (COVID-19) and contributes to disease severity has not been well characterized. In this study, we used bulk and single-cell RNA-sequencing data from the blood of severe COVID-19 patients to investigate the role of ERVs in the immunopathology of COVID-19. Upon quantification of proviral ERV expression, we identified 33 independently transcribed ERV loci that were differentially expressed in severe COVID-19 compared with healthy individuals (severe COVID-19 signature ERVs). ERV activation was associated with changes in the epigenetic regulators of ERVs and strongly correlated with key inflammatory pathways of COVID-19, including neutrophil degranulation, interleukin signaling, and inflammasome activation. Six of the upregulated ERV loci were specific to intensive care unit (ICU) admission and significantly correlated with disease severity, as measured by hospital-free days at day 45 (HFD-45). Finally, at the single-cell level, we detected significant upregulation of severe COVID-19 signature ERVs in erythroid-like and erythroid precursor cells and macrophages of patients with severe disease. Within these cell populations, we found evidence of ERV-associated differences in inflammatory gene expression, whereby cells expressing signature ERVs showed heightened expression of innate immune genes compared with cells not expressing these ERVs. Together, our study unmasked specific ERV loci activated in severe COVID-19 that are linked to the immunopathology that drives progression of severe COVID-19.

15
Characterization of the IGH locus and tissue specific immunoglobulin repertoires in turbot (Scophthalmus maximus).

TOUCEDO, R.; Zhu, Y.; Moledo, S.; Gambon Deza, F.; Boudinot, P.; Santos, Y.; MAGADAN, S.

2026-07-03 immunology 10.64898/2026.06.30.735498 medRxiv
Top 0.2%
1.0%
Show abstract

Turbot (Scophthalmus maximus) is an important aquaculture species, but the genomic organization and expressed diversity of its antibody repertoire remain incompletely characterized. In this study, we annotated the immunoglobulin heavy chain (IGH) locus using the haplotype resolved fScoMax1.1 genome assembly, and we used this as a reference to profile the expressed turbot IgM, IgD and IgT repertoires in skin and spleen. The primary IGH locus was located on chromosome 19, spanned approximately 72 kb, and contained 25 IGHV genes, including 24 functional genes and one pseudogene, together with three IGHD, seven IGHJ and three IGHC genes corresponding to IgT, IgM and IgD. Comparison with the alternate fScoMax1.1 haplotype and a second turbot genome assembly showed conserved IGHD, IGHJ and IGHC content, whereas IGHV gene number differed among assemblies. High throughput 5RACE repertoire sequencing revealed isotype and tissue associated differences in expressed IGH diversity. IgM represented the dominant productive repertoire in both skin and spleen and showed the highest clonotypic diversity, particularly in spleen. IgD displayed an intermediate profile, whereas IgT was more enriched in skin and exhibited the strongest clonal restriction. IGHV subgroup usage was dominated by IGHV3 in IgM and IgD, whereas IgT showed a distinct profile characterized by preferential use of IGHV4, especially in skin. Gene level analysis further showed broad IGHV-IGHJ pairing in IgM and IgD, with preferential use IGHJ3 segment, while IgT sequences paired exclusively with IGHJT. Clonotype sharing between skin and spleen was isotype dependent, being strongest for IgT, intermediate for IgM, and negligible for IgD, suggesting that clonal expansion did not necessarily predict inter tissue trafficking. Together, these results provide a curated genomic and expressed repertoire framework for turbot IGH genes and reveal isotype specific organization of antibody diversity, with IgT displaying a particular repertoire pattern.

16
Multiomics profiling of APECED peripheral blood highlights compositional increases in alternatively activated B cell subsets

Galbavy, W.; Kim, H.; Klotz, B.; Malbec, M.; Tasker, C.; Devlin, J. C.; Daniel, B.; Lim, W. K.; Benitez, A. A.; Haxhinasto, S.

2026-07-27 immunology 10.64898/2026.07.24.740518 medRxiv
Top 0.2%
1.0%
Show abstract

APECED (Autoimmune PolyEndocrinopathy Candidiasis Ectodermal Dystrophy) is a rare syndrome of multi-organ autoimmunity driven by the presence of self-reactive T cells and autoantibodies caused by mutations in the gene Autoimmune regulator (AIRE). Compared to the well-defined role of AIRE in establishing and maintaining T cell central tolerance, less is known about how AIRE deficiency impacts B cell phenotypes that may contribute to a breakdown in peripheral B cell tolerance. Here we analyzed serum and peripheral blood cells from APECED patients and healthy donors using autoantibody profiling, proteomics, flow cytometry, scRNAseq based subset analysis, BCRseq, and autoantigen binding assays finding significant changes to the APECED B cell compartment. We show that while Naive and Transitional B cells are reduced, alternatively activated B cell subsets are expanded in APECED patients including IgM CD27+ and class switched Atypical B cells which exhibit BCR chain features prone to autoreactivity. Antibodies derived from either APECED or Healthy class switched atypical B cells bind autoantigens at significantly higher rates than control IgG B cells. Serum autoantibodies and proteomics highlight APECED common and patient variable changes. These results together show that APECED causes a compositional shift towards subsets that may promote broken B cell tolerance and autoimmunity.

17
Single-Cell Transcriptomic Analysis of the Immune Response to CHIKVInfection

Huang, P.; Wang, H.; Xu, S.; Li, M.; Guo, M.; Wang, H.; Gou, X.; Wang, C.; He, Y.; Pan, W.

2026-06-30 immunology 10.64898/2026.06.24.734421 medRxiv
Top 0.2%
1.0%
Show abstract

Abstract Background: The Chikungunya virus (CHIKV), a re-emerging mosquito-borne alphavirus, is responsible for acute febrile illness and severe polyarthralgia. Although both innate and adaptive immune responses influence the disease outcomes, the detailed cellular immunopathogenesis of CHIKV in peripheral blood is not yet fully elucidated. Methods: We conducted single-cell RNA sequencing (scRNA-seq) on peripheral blood mononuclear cells (PBMCs) obtained from patients acutely infected with CHIKV and from healthy control subjects. Cellular interactions were inferred, and the transcriptomic results were orthogonally validated through quantitative real-time PCR (qPCR) and Enzyme-Linked Immunosorbent Assay (ELISA) to assess systemic interferon-stimulated responses. Furthermore, a comparative analysis was performed using publicly available single-cell data from Dengue virus (DENV) infections. Results: CHIKV infection significantly altered the immune system, increasing monocytes and dendritic cells while reducing T and B lymphocytes. Monocytes and NK cells showed strong activation of interferon-stimulated genes (ISGs). Monocytes were identified as key in driving inflammatory and immune responses. In adaptive immunity, CHIKV led B cells to become plasmablasts with antiviral immunoglobulins and caused T cells and NK-like T cells to show signs of cytotoxicity and exhaustion. Validation showed increased levels of IFN-{gamma}, IFN-{beta}1, MX1, and ISG15. CHIKV triggered a more intense, monocyte-driven interferon response than DENV. Conclusions: Acute CHIKV infection induces a systemic interferon response predominantly centered on monocytes, accompanied by significant alterations in adaptive immunity. Circulating ISG products, including MX1 and ISG15, reflect the transcriptomic activation and may serve as potential biomarkers for assessing the early intensity of innate antiviral responses.

18
Is APOE ε2 always a protective allele? Deviations in Hardy-Weinberg equilibrium in admixed Brazilian elderly individuals

Santos, G. d. N.; Rodrigues, P. H. S.; Passos, C. H.; Paco, S. L. G.; Ignacio, I. B.; de Alexandria, M. A. L. S.; Bastos, A. O.; Veronezz, L. A.; Neto, F. A. d. O.; Bardella, M. U.; Suemoto, C. K.; Leite, R.; Meyer, D.; Grinberg, L.; Naslavsky, M. S.

2026-06-23 genetics 10.64898/2026.06.20.733520 medRxiv
Top 0.2%
0.9%
Show abstract

The APOE gene is a critical determinant of human healthspan and longevity, with the rare{varepsilon} 2 allele traditionally viewed as a universal protective factor against Alzheimer s disease (AD) and a driver of exceptional lifespan. However, this protective paradigm is predominantly derived from European-centric cohorts, leaving the evolutionary and clinical impacts of{varepsilon} 2 across diverse, highly admixed populations largely unknown due to a lack of local ancestry (LA) resolution. To investigate how local genomic backgrounds modulate APOE survival dynamics we analyzed two Brazilian sample collection of older adults from Sao Paulo city: the Biobank for Aging Studies (BAS, n = 716), a post-mortem autopsy study of naturally deceased individuals; and the Health, Well-being and Aging Study (SABE, n = 952), a census-based elderly sample collection. We evaluated deviations from Hardy-Weinberg equilibrium (HWE) using robust permutation-based models to capture ongoing selective and mortality pressures at the APOE locus. While global APOE frequencies adhered to HWE, integrating LA unveiled striking, mirrored ancestral deviations. Our findings reveal that APOE {varepsilon}2 homozygotes with African ancestry significantly contribute to deviations from HWE in the BAS, with an excess of {varepsilon}2AFR/{varepsilon}2AFR homozygotes observed (p = 0.0196). These distinct HWE deviations demonstrate that an African LA background acts as a genetic buffer, attenuating the phenotypic extreme effects of APOE alleles. Furthermore, we observed an excess of the{varepsilon} 4 European haplotypes in the BAS, which is consistent with a mortality pressure allelic effect in the European LA context. Conversely, the{varepsilon} 4AFR/{varepsilon}4AFR combination was overrepresented in the SABE. While this buffering mechanism mitigates{varepsilon} 4 toxicity, it simultaneously dampens the exceptional longevity advantage typically conferred by the{varepsilon} 2 allele, leading to its neutral accumulation in the post-mortem cohort. Our study challenges the "one-size-fits-all" assumption of APOE biomarkers, demonstrating that{varepsilon} 2 protective mechanisms are context-dependent and modulated by local genomic backgrounds in admixed populations.

19
Characterisation of Peripheral Blood B Cell Receptor Repertoire in Severe Eosinophilic Asthma and EGPA

Arora, J. K.; Bessell, E.; Beyatli, S.; Thenet, D.; Brown, J.; Nissim, A.; Lewis, M. J.; James, L. K.; Pfeffer, P. E.

2026-06-20 immunology 10.64898/2026.06.16.732558 medRxiv
Top 0.2%
0.8%
Show abstract

BackgroundSevere eosinophilic asthma (SEA), eosinophilic granulomatosis with polyangiitis (EGPA) and nasal polyposis (NP) are immune-mediated diseases characterised by eosinophilic inflammation. However, there is also increasing interest in the potential pathological roles of autoantibodies in these diseases. Understanding their B cell receptor (BCR) repertoires may provide valuable insights into disease mechanisms, and potential role of B cells in their pathology. MethodsWe conducted BCR repertoire sequencing using peripheral blood from 43 patients, comprising SEA with nasal polyps (SEA+NP), SEA without nasal polyps (SEA-NP), and EGPA, along with 16 healthy controls (HCs). ResultsCompared to HCs, patients with EGPA exhibited increased relative proportions of IgA1, IgG1, IgG2, and IgG4 subclasses. Similarly, SEA-NP patients demonstrated significantly high proportion of IgG2 sequences. Notably, the IgG4 subclass was significantly elevated across all patient groups compared to HCs. Patients receiving anti-IL-5/5R biologic treatments showed increased relative proportions of IgA2 and IgG2 subclasses compared to untreated patients. Some variation across participant groups in mean somatic hypermutation and mutation frequency was evident. 1,508 clones shared across patients, but not healthy controls, were evident though the majority showed low clonal expansion. Nevertheless, a few shared clones did show either high prevalence across patients and/or higher clonal expansion. ConclusionChanges in BCR repertoires in SEA/EGPA are consistent with a pattern of a more mature B cell component in the periphery and with the T2 inflammatory response observed in SEA and EGPA. BCR clonotypes shared across patients were evident, however, whether such clonotypes are pathological in SEA/EGPA requires further investigation.

20
Global population frequencies of NAT2 star alleles observed in three large biobanks

Sangkuhl, K.; Whirl-Carrillo, M.; Woon, M.; Venkatesh, R.; Keat, K.; Whaley, R.; Ritchie, M. D.; Klein, T. E.

2026-06-11 genetic and genomic medicine 10.64898/2026.06.09.26355281 medRxiv
Top 0.2%
0.8%
Show abstract

NAT2 is an important pharmacogene which encodes the N-acetyltransferase 2 enzyme that is involved in the metabolism of multiple medications, and variants in this gene can affect patient response to these medications. CPIC has published a clinical guideline for prescribing hydralazine using NAT2 genotypes. Just prior to the guideline, updated NAT2 star allele numbering and definitions were released, differing somewhat from the historical nomenclature. Clinical pharmacogenomic testing panels often test for the most common star alleles, so knowledge of the most common updated NAT2 star alleles is critical for the implementation of the CPIC NAT2/hydralazine guideline. We first determine NAT2 diplotype frequencies from UK Biobank (UKBB) 200k phased genomes, then analyzed allele, diplotype, and phenotype population frequencies from the All of Us Research program, PennMedicine BioBank (PMBB) and UKBB 500k datasets. We found that analyzing NAT2 diplotypes from phased data provides critical information for algorithms designed to predict diplotypes from unphased data. We observed that NAT2*5, *6, and *4 were the most common star alleles in that order, and the top 11 most frequent NAT2 star alleles were the same across all biobanks. However, differences in star allele frequencies across biogeographical populations were observed. The largest difference led to a higher frequency of NAT2 poor metabolizer phenotypes as compared to rapid and intermediate metabolizer phenotypes in all global populations except in the EAS population, where NAT2 poor metabolizers were in the minority.